Relation Between Bulk and Interface Descriptions of Alloy Solidification 



Alexander L. Korzhenevskii 
Institute for Problems of Mechanical Engineering, RAS, 
Bol'shoi prosp. V. 0., 61, St Petersburg, 199178, Russia 

Richard Bausch 

Institut fur Theoretische Physik IV, Heinrich-Heine-Universitat Diisseldorf, 
Universitatsstrasse 1, D-40225 Diisseldorf, Germany 



(N 

o 

(N 



CZ2 



Rudi Schmitz 

Institut fur Theoretische Physik C, RWTH Aachen University, Templergraben 55, D-52056 Aachen, Germany 

(Dated: July 3, 2012) 

From a simple bulk model for the one-dimensional steady-state solidification of a dilute binary 
alloy we derive an interface description, allowing arbitrary values of the growth velocity. Our 
derivation leads to exact expressions for the fluxes and forces at the interface and for the set of 
Onsager coefficients. We, moreover, discover a continuous symmetry, which appears in the low- 
velocity regime, and there deletes the Onsager sign and symmetry properties. An example is the 
generation of the sometimes negative friction coefficient in the crystallization flux-force relation. 
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INTRODUCTION 

The interface kinetics at the solidification front of a 
growing dilute binary alloy has been described by Baker, 
and Cahn [l[ within the framework of linear irreversible 
thermodynamics. An apparent violation of the Onsager 
symmetries in their approach has led to several partly 
controversial discussions in papers by Caroli, Caroli, and 
Roulet 1 2], by Kaplan, Aziz, and Gray |3|], and in the 
review y| by Hillert. Another apparent paradox is the 
possible appearance of a negative friction coefficient in 
the crystallization flux-force relation. This unexpected 
effect has been observed in the phase-field approach by 
Karma and Rappel [H, and in Ref. |6| has been shown 
to be compatible with the principle of a positive entropy 
production. A thorough discussion of this point is also 
contained in a recent phenomenological approach to the 
one-dimensional steady-state solidification by Brener and 
Temkin [7|. 

In the present work we derive the interface description 
for the solidification of a dilute binary alloy from a simple 
bulk model. The latter is a conveniently chosen specific 
version of the capillary-wave model, derived in Ref. [6| 
from a phase-field model. As in Ref. Q, we consider 
the scenario of a one-dimensional steady-state motion of 
the solidification front, however, without any restriction 
on the value of the growth velocity. Our main results 
include exact expressions for the fluxes and driving forces 
at the interface, carrying the full velocity dependence 
due to the solute trapping effect. We also obtain exact 
expressions for the Onsager coefficients, which, in the 
basis of crystallization and diffusion fluxes and forces, 
turn out to be independent of the growth velocity. 

Another result of our investigations is the observation 



that a continuous symmetry arises in the low-velocity 
regime due to a linear dependence between the two fluxes 
of the model. This symmetry conserves the transport 
equations through the interface and, as a consequence, 
the entropy production within the interface region, but 
generates two independent redundant parameters in the 
set of transport coefficients. In general, this effect deletes 
the Onsager sign and symmetry conditions of this set, 
and, by properly choosing the parameters, leads to the 
appearance of a sometimes negative friction coefficient in 
the crystallization flux-force relation. As to be expected, 
this apparent anomaly can easily be seen not to violate 
the principle of a positive entropy production. 



CAPILLARY- WAVE MODEL 

Based on a phase-field description, we have established 
in Ref. @ a capillary-wave model for the solidification of 
a dilute binary alloy where the interface position Z(x, t) 
enters as an extra field variable in addition to the solute 
concentration C(r, i). In the one-dimensional case the 
effective Hamiltonian of our model reads 
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dz[C(z,t) - U(z - Z{t))f 



(1) 



Here, in addition to the coupling constant k, the function 
XJ{z — Z) enters as an input quantity, which, following 
from the condition SH/SC = 0, means the equilibrium 
concentration of the solute. The equations of motion of 
the model are given by 
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d t C = d z D{z - Z(t))d, 
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where A and D(z — Z) are kinetic coefficients, measuring 
the mobility of the interface and the local diffusivity of 
the solute. Instead of including an external driving force 
in the first of these equations, as done in Ref. [hL the 
solidification process will be activated, as in Ref. [7j, by 
imposing appropriate boundary conditions on the second 
of the equations ((2|). 

In a steady-state solidification process with constant 
growth velocity V in z-direction, the solute concentration 
in the comoving frame only depends on the variable 



C = z - Vt 



(3) 



As a result, insertion of Eq. ((TJ) into the first equation in 
Eqs. © yields 



V = -Ak 



d(u'(0 [ C (0-u(0] 



(4) 



The second equation in Eqs. ([2]) can be integrated once, 
which, with the boundary condition 



C(-oo) = Co , (5) 
leads to the first-order differential equation 



[C -C(0}V = D(O[C'(O-U'(O] 



(6) 



Eq. ^ can be rewritten in form of the force-balance 
equation 



V (AC) 2 

-_ + F + G{v) + K L_L = (7) 

where the first term has the form of a viscous-friction 
force. The contribution 



AC = C L - C s 



(10) 



means the miscibility gap. Evidently, F describes the 
effect of a quench from the equilibrium value Cs to some 
non-equilibrium value Co, which for Cs — Co > induces 
a growth rate V > 0. The contribution 



G{V) = -k d(U'(OlC(()-C } 



(11) 



is a drag force, generated by the comoving solute layer at 
the interface. Using the steady-state boundary condition 



C(+oo) = C(-oo) , (12) 
it can be rewritten, due to Eq. ©, in the form 



G{V) = — k 



+°° n(f) 



(13) 



which, for velocities < V < oo, is negative and shows 
the appropriate behavior G(V) — > for V — > oo. At 
V = 0, corresponding to the choice C = Cg, Eq. (O 
implies G(0) = — k(AC) 2 /2, elucidating the role of the 
last term in Eq. (|7|). 

Ref. [6| also provides an expression for the density p(C) 
of the entropy production. The phase field entering 
this expression, is, according to Ref. 6], related to /7(C) 
by the equation 



*(0 - ^^(0 



Cl + Cs 
AC 



(14) 



and obeys the normalization condition 



^Ka c $(C)] 2 = i 



(15) 



where £ is a microscopic length scale of the order of the 
interface thickness. Due to Eq. (fl4")) , the related passage 
in the discussion of Ref. @ leads to the obviously positive 
expression 



^ = «(AC)(C S -C ) (8) 

acts as an external driving force where, in terms of the 
solute concentrations 

C s = U(-oo) , C L = U(+oo) (9) 
in the liquid and solid phases, 



1 V 2 ( 2 \ 2 1 2 

KC) = fi y -^MOf 

1D(C) ( 5H V 
+ T—{ d <5C(0) ■ (16) 

for the density of the entropy production at some fixed 
temperature T. 

In order to explicitly evaluate Eqs. ©, ©, and (fTl)l) . 
we now consider the specific model, defined by 
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and the identity 



U(0 = C S + -(C L -C S ) 


x < C + 2£ 
4£ 





4£ 



2£ < C < +2£ 
+ 2£ < C, 



£'(C)=£' s + -(Z> i -£» s ) 



- 2^ < c < +2e 

+ < C ■ 



(17) 



(18) 



Here, £>s and Dl are the diffusion coefficients of the 
solute in the solid and liquid phases, and the width 4£ of 
the well-encompassed interface region has been chosen to 
ensure the normalization (|15[) . In the special case Ds = 
the solution of Eq. ©J subject to the boundary condition 
(fl"2|) . has the simple form 



C-=C 



(22) 



implied by Eq. (fl~9|) . By insertion of the solution (fT9)) into 
the force balance ([7]) one, furthermore, finds the relation 



C = C s - (AC) 
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(23) 



V c V D + 2V, 
written in terms of crystallization and diffusion velocities 

V C = «(AC) 2 A , Vjj = ~7£P • (24) 

According to the result (|2"3"|) . the limit F — >• of the 
driving force © not only is compatible with the behavior 
V -> 0, but also with the possibility V ->■ (Vfc - Vb)/2. 
In the process of establishing the interface description, 
this observation leads us to take care of the full velocity 
dependence, appearing in the solution (|19l) . 

INTERFACE DESCRIPTION 



C(C) = Co + ^(C L - C S ) (l + 4^ ) ( I')! 



x < 





(C + 20 
4£exp 



c<-2e 
2^ < c < +2e 

+ 2$ < c, 



showing, due to the factor (1 + 4£V/Dl) _1 , the solute- 
trapping effect. Some properties of this solution have 
already been discussed in Ref. @ , and below will be used 
to establish the interface description of the solidification 
process. 

The entropy production P, generated in the interface 
region per unit area, follows from Eq. (|16j) by integration 
over the interface region — 2£ < ( < +2£. Taking into 
account Eqs. (TJ), ©, ©, (jTTJ) , and (|T5]). one obtains 



P = — 
T 

V 

+ T 

V 

T 



F + G(V) + -(AC) 2 



(20) 
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2£ 



dC[c - c(C)]«[c"(0 - u'(0\ 



F 



|(AC) 2 -|(C + -C_) 2 



where we have used the notation 



C± = C(C = ±20 , 



(21) 



The last expression in Eq. (|20| can be written in the 
standard form 



P — Ic'Ec + IdEd 



(25) 



where the crystallization and diffusion fluxes Ic,Id and 
the related forces Ec, Ed are exactly given by 



Ic = pV , 

I D = ±(C_-C + -AC)V , 



(C+ - CL - AC) , 



(26) 



T 



p meaning the density of the material. An equivalent but 
more familiar expression for Ed is 



E D = 5 ± ee 1^+20 -K-2£)] (27) 
where /i(£) is the chemical potential, defined by 



SH 



k[C(0 - U(Q] 



(28) 



The constitutive equations, connecting the fluxes and 
forces (|26l) . are obtained by evaluating Eqs. ^ and ([7]) 
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by means of Eqs. ([17]). (|T5| . and (HU). This leads to the 
result 



E c 
E D 



Rcc Rcd 
Rdc Rdd 



Ic 
Id 



(29) 



with the again exact expressions of the matrix elements 



V P 

AC AC 

Rcd — 2 , Rdc — 2 

P P 



(30) 



We notice that the fluxes and forces (|26p carry the 
full velocity dependence, reflecting the solute trapping 
effect whereas the kinetic coefficients (|3T)|) turn out to be 
independent of the velocity. The latter, moreover, obey 
the basic properties 



Rcd — Rdc > Rcc > , Rdd > , 
Rcc Rdd — Rcd — 



(31) 



of general Onsager coefficients. 

In principle, the results (|26|) . (|29|) . (|30]) are of limited 
interest, since the solution (|I 9|) . in combination with the 
resulting explicit form of the force balance (O, describe 
the complete steady-state behavior of our model in the 
interface and bulk regions. However, these results are a 
good starting point for an exhaustive discussion of the 
interface description in the often considered low- velocity 
regime. 

LOW- VELOCITY LIMIT 



with coefficients 



Sec 



S, 



CD 



s 



DC 



Sdd = M 

AC 
P > 



(35) 



leads, due to Eq. (|3"2"|) . to separate cancellations of the 
contributions containing the free parameters A and /i in 
Eq. (|29p. However, the matrix elements ([3T))) transform 
into the A, /i-dependent coefficients 



Iter = I ^ + 



Rcd 



V c 

(A + 2)^ 

P 



(f 

Rdc = 



Rdd = p + 4 , 
(36) 



P 



We mention that this result is in accordance with an old 
observation by Hooyman and De Groot [8|. It should 
also be noted that in the subdomain 



A = (j, > 



(37) 



the coefficients (|3^)l obey the Onsager conditions (|3"Tj) in 
a formal sense. In this domain they depend, however, on 
the redundant parameter A. 

Generally, it is fairly allowed to ignore the constraint 
(1371) and to choose in Eqs. (|36l the values 

A = /i = -2 , (38) 

which leads to a diagonal form of the matrix R and, 
consequently, to the decoupled transport equations 



Whereas the fluxes, defined in Eqs. (|26l) . are linearly 
independent for arbitrary values of V, one finds in linear 
order in V the connection 



Id — ic ■ 

P 



(32) 



This relation implies the existence of a two-parameter 
continuous symmetry, which changes the coefficients 
([3T))) . but conserves the kinetic equations (|2"9"|) . and, as 
a consequence, also the entropy production (|25[) . In fact, 
the transformation 



R = R + S 



(33) 



involving the matrices 



K 2 , 



Er 



id 



(39) 



explicitly showing that the coefficient, connecting Ec and 
Ic can become negative. Due to Eqs. (|25l) and (|3"2"j) the 
related entropy production is, in the presently discussed 
low-velocity regime, given by the expression 



1 (—] ( — — 
T\~) \Vc + Vd~ 



l c > 



(40) 



which obviously is positive. 



TWO-COMPONENT REPRESENTATION 



R 



Rcc Rcd 
Rdc Rdd 



S = 



S q cc S q CD ) , (34) 

ODC JDD ' 



The fluxes, associated with the two components of the 
solidifying alloy, are given by the relations 
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(£)-(! .- x x)(f D ). <-) 

where, by definition, 

X = ^-(C- + C+ + AC) . (43) 
2p 

In the entropy production 



The last statement seems to contradict a conclusion 
by Baker and Cahn which later has been criticized 
by Caroli, Caroli, and Roulet @ as a fallacy, resulting 
from the assumed absence of solute diffusion in the solid 
phase. Contrary to this claim, the erroneous conclusion 
in Ref. [l[ is, to our mind, due to an algebraic mistake, 
which we here clarify, picking up from Ref. [l[ a model, 
which has only diagonal elements in the inverse relation 
ofEq. (47), 

J A = L AA F A , J B = L BB F B . (49) 

From Eqs. (|H) . (|42|) . (|49) one correctly obtains the first 
constitutive equation in the C, D basis, 



P = F a J a +FbJb 



(44) 



arising from the similar expression (|2"51 , all X-dependent 
terms cancel. The special choice (|4"3") . however, leads to 
the relations 



J A = [p-C-)V , Jb = C-V 



(45) 



where, parallel to Eq. ([28]) . 



MO = - 
p 



i - c(0 



ac(0 



-ic(0-u(()} 2 



(46) 



which, together with Eq. (28), provides agreement with 
standard notations. It is obvious that, in discussions of 
the low- velocity limit, one would replace the quantity 
(|43| by the approximate value X = (AC) /p. 

which together with Eq. (|2"51) establishes agreement 
with the notions, used in Ref. I]. It is obvious that, in a 
discussion of the low- velocity limit, one would replace the 
quantity (|4"3"| by the approximate value X — (AC) /p. 

The constitutive equations (|29l) , finally, transform into 



F b 



Qaa Qab 
TVd V Qbc Qbb 



Ja 
Jb 



(47) 



where the kinetic coefficients in this basis are related to 
those in the C, D basis by 



Qaa Qab 
Qbc Qbb 
1 -X 
I I -X 



(48) 



Rcc Rcd 
Rdc Rdd 



I 1 

-X 1 - X 



One easily verifies that the Onsager symmetry and sign 
conditions (|3"Tj) also are valid in the A, B basis. 



(Laa + Lbb)Ec 

+ [(1-X)L B b-XLaa]E d 



(50) 



In order to attain the second equation, which determines 
Id, the authors of Ref. [l[ just use the relation (|3"2")l 
together with the result (|50p . This does not exhaust the 
full content of Eqs. (f49|) . and as a result, their mapping 
from the A, B to the C, D basis is singular instead of 
being a one-to-one mapping. The correct procedure is, 
to calculate from Eqs. (|4"Tj) the flux 



I D = (1 - X ) J B - XJ A , 
which, together with Eqs. (|4TJ)) . yields 



I D = [(l-X)L 



BB 



XLaa]Ec 
\{1-X) 2 Lbb + X 2 L A a]E d 



(51) 



(52) 



Comparing this with Eq. (|50p , one immediately observes 
that the diagonal matrix elements are positive and that 
the off-diagonal matrix elements are symmetric. It is a 
simple exercise, to verify that the determinant, formed 
with these elements, also is positive, in accordance with 
our general statement below Eq. (|4"5)l . 



CONCLUSIONS 

The model, leading to the exact results (26), (29l), (30]) 
presumably is the simplest one, which can exactly be 
evaluated. We also have considered the case of a nonzero 
diffusion constant D$ in the model equation (|I8[) , which 
only led to slightly more complicated expressions for the 
Onsager coefficients. 

One of the main insights, drawn from our analysis, 
concerns the redundancies in the transport coefficients, 
which appear in the low-velocity regime, if one tries to 
sustain the formal structure of the Onsager approach. 
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Another conclusion is that, in the discussion of steady 
states, one should not overlook the possibility of a finite 
growth velocity, which is clearly compatible with a small 
driving force. This case is necessarily connected with 
the appearance of a negative friction coefficient in the 
crystallization flux-force relation, and completely covered 
by our analysis. 
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